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The invention relates to a liquid crystal based polarimetric system 
for analyzing a sample, a process for the calibration of this polarimetric 
system and a polarimetric measurement process. 

In order to measure parameters which are representative of a 
5 sample (for example, of its composition and thickness), it is conventional 
to make use of an ellipsometer. 

Ellipsometry is a powerful non invasive probe in which reflectance 
or transmittance data are measured by electromagnetic radiation 
outcoming from a sample. Briefly, the surface of a studied sample is 
10 illuminated by a photon beam that is reflected or transmitted and the 
polarization state of the outcoming beam is compared to that of the 
incident beam. 

This conventional ellipsometry method proves satisfactory when 
the reflected light is still totally polarized (even though elliptically), which 
15 is indeed the case, among other examples, for isotropic layers with 
smooth interfaces. Such samples, which can be considered as “dichroTc 
retarders” (DR) are usually characterized by the ellipsometric angles 
(T,A) defined by 

20 rp/rs= tan(T) exp(iA) (1) 

where rp and rs are respectively the amplitude reflectivities of the sample 
for linear polarizations in the incidence plane (rp), or perpendicular to this 
plane (rs). 

25 However, to study samples that cannot be described as DRs, such 

as partially depolarizing materials, a more general method is required. 

Polarimetric systems enable to measure all the polarization 
components of light in any sample. 

The polarization state of light can be represented in the most 
30 general case by a four-dimensional vector, called the Stokes vector (S). 

A description of this can be found in the work of Azzam and 
Bashara entitled “Ellipsometry and polarized light”, North-Holland, pp. 55- 
60. 

The Stokes vector consists of the I, Q, U and V coordinates. I 
35 stands for the total intensity, while the other components are defined as 
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the differences between the intensities measured through different pairs 
of “complementary” polarizers (linear polarizers set vertical and horizontal 
for Q, at + 45° and - 45° azimuthal angles for U, and left and right 
circular polarizers for V). 

5 The interaction of light with any sample can then be represented 

by a matrix, so-called, the Mueller matrix, of dimensions 4x4 with 
therefore 16 real coefficients. 

The extraction of the 16 parameters during polarimetric 
measurements provides a complete characterization of the said medium. 
10 For a DR characterized by ellipsometric angles ('T, A) (see eq. (1)) 

the Mueller matrix elements are the following: 

- upper diagonal block 

Mi-i=M22“'^i Mi2“M2i~“ t cos(2'P) 

- lower diagonal block 

15 M 33 =M 44 =x sin(2T) cos(A) , M 34 =-M 43 = X sin(2T) cos(A) 

- other elements : My = 0 (2) 

where the additional parameter x is proportional to the overall intensity 
transmission or reflection coefficient of the sample. We point out that this 
Mueller matrix has 

20 - two real eigenvalues : (Ari=2x sin^(40, ^=2x cos^(TO), 

- and two complex conjugate eigenvalues: 

(2ci=x sin(2T0 exp(iA), Ac 2 =x sin(2^ exp(-iA)). (3) 

Many designs have been proposed, and demonstrated, for Mueller 
polarimetric systems. All of them comprise a polarization state generator 
25 (PSG) which modulates the Stokes vector (Sin) of the light impinging on 
the sample and a polarization state detector (PSD) which analyzes the 
polarization (Sout) of the light outcoming from the sample. It is customary 
to define the modulation matrix W as a 4x4 real matrix whose columns 
are the four Stokes vectors Sin generated by the PSG. Reciprocally, the 
30 four dimensional signal vector D eventually delivered by the PSD is 
related to the Stokes vector Sout of the light outcoming from the sample 
by a linear relationship D = A Sout. where A is the (4x4 real) analysis 
matrix representing the PSD. A raw measurement actually consists of 16 
values of the signal, which form a matrix B = A M W, where A and W are 
35 respectively the analysis and modulation matrices defined above, and M 
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the Mueller matrix of the sample. If A and W are known, M can be 
extracted from the raw data B as M = A'^ B W\ The determination of A 
and W is called calibration of the polarimeter. Clearly, the instrument 
must be designed in such a way that A and W are non singular. 

5 Moreover, in order to optimize error propagation from the raw 
measurement B to the final result M, the analysis and modulation 
matrices A and W must be “as close as possible” to unitary matrices. The 
best criterion in this respect is to optimize their condition numbors s(A) 
and s(W), which are the ratios of the smallest over the largest of their 
10 singular values {see for example Press W.H., Teukolsky S.A., Vetterling 
W.T. and Flannery B.P.; Numerical Recipes in Fortran, Cambridge 
University Press, p 53, who actually define the condition number as the 
reciprocal of that defined here; see also Compain E. and Drevillon B., 
Rev. Sci. Instrum. 69, (1998) 1574}. 

15 In a PSG, the light polarization can be modulated by a variety of 

devices such as discrete components inserted and then removed from 
the light path {Bickel W.S. et al.; Am. J. Phys 53 (1984) 468}, rotating 
retardation plates (Goldstein D.H.; Appl. Opt. 31 (1992) 6676}, rotating 
compensators (Collins R.W. and Koh J.; J. Opt. Soc. A 16, (1999) 1997}, 
20 Pockels cells (Delplancke F.; Appl. Opt. 36 (1997) 5388 and Compain E. 
and Drevillon B.; Rev. Sci. Instrum. 68 (1997) 2671} or photoacoustic 
modulators (Compain E. and Drevillon B.; Rev. Sci. Instrum. 69, (1998) 
1574}. For PSD, one can use the same devices and a single detector, or 
a “parallel” analysis of light polarization through polarization-sensitive 
25 beamsplitters and simultaneous measurement of the separated beams 
by several detectors (Azzam R.M.A., Opt. Acta 29 (1982) 685, 
Brudzewski K.; J. Modern Optics 38 (1991) 889, Khrishnan S.; J. Opt. 
Soc. Am A 9 (1992) 1615, Compain E. et al.. Thin Solid Films 313 

(1998)}. 

30 This variety of designs leads to a variety of characteristics, some 

of which are not compatible with each other; for example, high frequency 
polarization modulation by resonant devices such as photoelastic 
modulators allows efficient rejection of low-frequency noise, but for 
imaging applications by slow detectors such as typical Charge Coupled 
35 Devices (CCD), a stroboscopic illumination is then needed. Another 
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important requirement for imaging applications is that the polarimetric 
elements exhibit large enough acceptance angles together with small 
enough aberrations. As a result, the devices which best meet these 
requirements are those based on low order retarding plates, which are 
5 either rotated {Pezzaniti J.L. and Chipman R.A., SPIE proceedings 2297 
(1994) 468} or inserted and removed between sequential measurements, 
or those based on liquid crystal (LC) variable retarders. 

Liquid crystal cells (LC) are electrically controlled low order 
retardation plates. Two types of such devices are currently available. 
10 First, nematic liquid crystals (NLC) provide variable retardation with fixed 
orientation of slow and fast axes, with typical response times of the order 
of 10 to 100 ms. On the other hand ferroelectric liquid crystals (FLC) 
provide fixed retardation, but with slow and fast axis directions which can 
be electrically switched between two azimuthal angles separated by 45°, 
15 in times typically shorter than 100 ps. 

These elements provide non resonant polarization control which is 
naturally well suited for polarimetric imaging by a CCD. Therefore, 
devices comprising liquid crystal cells have been proposed for 
polarimetric imaging within the frame of conventional ellipsometry, i.e. for 
20 samples behaving as DR {Oldenbourg R. et al.; J. Microscopy 180 (1995) 
140} and led to commercially available devices (Pol-Scole, by CRI, Inc. 
Boston). 

Stokes polarimetry, i.e. polarimetry performed using a sole 
polarization state detector and no polarization state generator has also 
25 been performed, essentially in solar astronomy. The device consisted of 
two nematic LCs followed by a linear polarizer {Hofmann A.; SPIE 
proceedings 4133 (2000) 44} or even more complex assemblies, 
including for example two ferroelectric LCs, two fixed X/8 retardation 
plates and a linear polarizer {Gandorfer A.M.; Opt. Engineering 38 (1999) 
30 1402} or one ferroelectric, two nematic LCs and two quarter wave 

retardation plates {November L.J. and Wilkins L.M.; SPIE proceedings 
2265, 210}. 

An imaging Mueller polarimeter has been realized by using 
nematic LC cells {Bueno J.M.; J. Opt. A; Pure and Applied Optics 2 
35 (2000) 216}. In this device, the PSG and PSD have the same design: 
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each of them consists of one LC and one quarter-wave plate, the latter 
being mechanically inserted in and removed from the light path between 
acquisitions of raw images. This device has been used for polarimetric 
imaging of human eye, including retina and cornea. 

5 However, the Mueller polarimeters described above, suitable for 

polarimetric imaging by slow devices such as CCDs, present two 
significant shortcomings. 

First, their operation involves mechanical motion of optical 
elements, which are either rotated or moved in and out of the light path. 

10 Second, their calibration relies on the characterization of individual 

optical components (polarizers and retardation plates). As a result, the 
accuracy of the overall calibration of the polarimeter is limited by the 
accumulation of the errors on the knowledge of each of those 
components, and on their positions. Furthermore, the instrumental 
15 configurations are most frequently defined in such a way that the 
polarization states generated by the PSG, or those “filtered by the PSD 
are “simple” polarization states, such as linear (vertical and horizontal) 
and circular states, to reduce “cross-talk” and facilitate the overall 
calibration of the system. These “simple” configurations are far from 
20 those providing the highest values of s(A) and s(W), implying that for a 
given input noise on the raw data (B), the noise on the final result (M), is 

far from being optimized. 

A purpose of the present invention is to remedy the shortcomings 
mentioned above and to propose a polarimetric system having one or 
25 more of the following features and advantages : namely, polarization 
modulation by liquid crystals only, with no mechanically moving parts, 
providing a wide enough acceptance angle for typical imaging 
applications, a simple and compact design including optimization with 
respect to error propagation, and a fast data acquisition, to be usable for 
30 measurements in real time. 

Another purpose of the present invention is to provide two 
measurement processes, which can be used on the same instrument, 

and which yield: 

• In 3 simpUfiQd opQrstion mods: the classical ellipsometric angles 
('P, A) of a sample assumed to be a DR. In addition, the validity of 
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this assumption (which depends on the sample homogeneity, 
roughness...) is tested automatically with no extra measurement, 
while with usual ellipsometer said test requires a mechanical 
rotation of the instrument output arm, 

5 • In a complete operation mode: the complete Mueller matrix (M) 

of any sample under study, either in transmission or in reflection. 

The optimization of the PSG and PSD configurations with respect to error 
propagation implies that the states generated by the PSG (and those 
“filtered” by the PSD) are not “simple” ones, such as linear or circular, 
10 and thus usual calibration methods are not really adequate. The 
invention includes therefore an objective calibration process for each of 
the two types of measurement cited above, said calibration processes 
being at once accurate, rapid and easy to implement. 

To this end, the invention concerns a polarimetric system for 

15 analyzing a sample comprising 

- an excitation section emitting a light beam, said excitation section 
comprising a polarization state generator containing a polarizer 
linearly polarizing the incident light beam along a direction of 
polarization (i), 

20 - an analysis section comprising a polarization state detector 

containing an analyzer, and detection means, 

- a processing unit. 

According to the invention, 

• the polarization state generator (PSG) comprise a first and a 

25 second liquid crystal elements LCj (j=1, 2) having, for each LQ 

element of the PSG an extraordinary axis making an angle 0j with 
respect to the direction of polarization (i) and a retardation 5j 
between its ordinary and extraordinary axes. Said liquid crystal 
(LCs) elements are placed after the polarizer and are equivalent to 
30 electrically controlled retardation plates, 

• the polarization state detector (PSD) comprise a first and a 
second liquid crystal elements LCj 0 = 1 . 2) having, for each LCj 
element of the PSD an extraordinary axis making an angle 0’j with 
respect to the direction of polarization (i) and a retardation 

35 5’j between its ordinary and extraordinary axes. Said liquid crystal 
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LCj elements equivalent to electrically controlled retardation plates 
are positioned in reverse order in the PSD with respect to the LCj 

elements of the PSG. 

According to various embodiments, the present invention also concerns 
the characteristics below, considered individually or in all their technical 
possible combinations. 

• In a first embodiment, the polarization state generator (PSG) 
and the polarization state detector (PSD) comprise each a first and 
a second nematic liquid crystal elements NLCj (j= 1, 2). For each 
NLCj element of the PSG (respectively, for each NLCj element of 
the PSD), the extraordinary axis makes a fixed azimuthal angle 6j 
(resp. 0’j) with respect to the direction of the input polarizer of the 
PSG, (resp. the output analyzer of the PSD). The retardation 5j 
(resp. 6’]) between NLCi ordinary and extraordinary axes is varied 

by means of an electrical control, 

-The azimuthal angles 0’j are equal to 0j G=1, 2) and the 
retardations 8’j are equal to -8j G=1, 2) (modulo 2ti) for 
simultaneous optimization of the condition numbers s(W) and 
s(A) of the modulation and analysis matrices. 

- By means of proper driving voltages, the couple of 
retardations (8i, 82) takes sequentially the following values : 
(Ai,Ai), (A,,A2), (A2 ,Ai), (A 2,A2), Where Ai and A2 verify the 
formulae (315“ + p 90°) and (135° + p 90“) respectively, where 

p is the same integer in both formulae, and 

- The angles 0i and 02 verify the formulae (s 27“ + q 90“) and (e 
72“ + r 90“) respectively where s = + 1 has the same value in 
both equations while q and r are any integer, with tolerances on 
the angles 0i and on the retardations Ai equal to +-10“ and +- 
20° respectively. With such tolerances, the condition numbers 
s(A) and s(W) are then always between the maximum value 

(equal to -^»0,58) and 0,3 which implies that the noise in the 

v3 

final matrix M (which is inversely proportional to s(A) and s(W) 
for a given noise on the raw data B), never exceeds twice its 

minimum value. 
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- For spectroscopic applications, i.e. for operation at a variable 
wavelength, a monochromator is piaced before the polarizer of 
the PSG, or after the anaiyzer of the PSD, and the values of 
the retardations (5i, 62 ) are kept within the boundaries specified 
5 above by simply tuning the amplitudes of the control voltages 

according to the waveiength passing through the 
monochromator. The currently available NLCs can be used 
from 400 nm to 1500 nm typically. 

• In an alternative embodiment, the polarization state generator 
10 (PSG) and the polarization state detector (PSD) comprise each a 

first and a second ferroelectric liquid crystai eiemehts FLCj Q= 
1,2). For each FLCj eiement of the PSG (respectively, for each 
FLCj element of the PSD), the retardation 6j (resp. 5’j) between 
FLCj ordinary and extraordinary axes is now fixed. For a given set 
15 of driving voltages the extraordinary axes of the FLC make a 

couple of azimuthal angles ( 0 i, 02 ) (resp. ( 0 ’i, 0 ’ 2 ) ) with respect to 
the direction of the input polarizer of the PSG (resp. the output 
analyzer of the PSD), and then these angles are set sequentially 
to (01,02), (0i+45°,02), (01,02+45°), (0i+45°,02+45°), 

20 -The azimuthal angles 0’j are equal to 0j 0=1, 2) and the 

retardations 5’j are equal to -5j 0 = 1 . 2) (modulo 2n) for 
simultaneous optimization of the condition numbers s(W) and 
s(A) of the modulation and analysis matrices. 

-The retardations (5i, 82 ) are given by 5i=80°+-15° and 
25 82 = 160 °+- 15 °, while the orientation angles (0i, 02 ) are given by 0i 

= 67°+-10° and 02 = 160°+-40°. With these values and 
tolerances, the condition numbers s(A) and s(W) are again 

between the maximum value (equal to -|=w0,58) and 0,3. 

-v/3 

- For spectroscopic applications, i.e. for operation at variable 
30 wavelengths, as the values of the retardations (5i, 82 ) are not 

electrically controllable as for the nematic crystal elements, the 
condition numbers s(W) and s(A) cannot be kept above 0.3 
throughout the visible with a PSG (or a PSD) comprising a linear 
polarizer and two FLCs only. However, with typical birefringence 
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dispersion of ferroelectric liquid crystais such a broadband 
optimization of the condition numbers (s(W) and s(A) > 0.3 in the 
whole spectrum covered by currently available FLCs (420 -800 
nm typically) is achieved by adding to the system another 
5 birefringent fixed element (retardation plate). An example of such 

an optimization is described hereafter, with a quartz plate 
inserted between the two FLCs. The broadband ferroelectric- 
based PSG then comprises: 

- a linear poiarizer, set at an orientation angle 0=0, 

10 - a first ferroelectric iiquid crystal, with a retardation 5i=90°+-5° at 

510 nm, set at an orientation angle 0i=-1O°+/-5“, 

- a quartz plate, providing a retardation 5q=90°+/-5° at 633 nm, 
set at an orientation angle 0q=5°+/-5°, 

- a second ferroelectric liquid crystal, with a retardation 

15 52=180°+/-15° at 510 nm, set at an orientation angle 02=71°+/- 

10 °, 

This polarimetric system is most conveniently coupled with a 
spectrometer placed after the analyzer of the PSD, and equipped 
with a multipoint detector (typically a CCD), which allows the 
20 polarimetric analysis to be carried out simultaneously in the 

whole spectral range defined by the currently available FLC 
elements, and which might be extended in the future with new 
FLC materials. 

Both embodiments described up to now are provided for illustrative 
25 purposes only and should not be used to unduly limit the scope of the 
present invention. For example, the polarimetric system is not limited to 
the use of NLCs or FLCs in the PSG and the PSD but a variety of 
devices combining ferroelectric and nematic liquid crystals can also 
be designed for simultaneous optimization of the condition numbers of 
30 the PSG and the PSD. 

• Said polarimetric systems are ellipsometers, 

• Said polarimetric systems are Mueller polarimetric systems for 
analyzing a sample represented by the sixteen coefficients of a 
Mueller matrix. 
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• The light beam emitted by the excitation section is in the 
spectral range 400-1500 nm for nematic liquid crystals and 420- 
800 nm for ferroelectric liquid crystals currently available, 

This spectral range might be extended in the UV or further into the 
5 IR with new LC materials with the different embodiments 

described within the scope of the present invention, 

• The excitation section comprises a monochromator positioned 
before the polarizer of the PSG, 

• The detection means comprises either a single detector, or a 

10 multipoint photosensitive detector, adapted with the processing 

unit to polarimetric imaging, 

• The multipoint photosensitive detector is a charge coupled 
detector (CCD), 

• For spectroscopic applications, the detection means may 

15 comprise a spectrometer, placed after the analyzer of the PSD, 

and preferably coupled with a CCD, to achieve polarimetric 
analysis simultaneously over the entire spectral range. 

The device can be used both in transmission and in reflection modes. 

The invention also concerns a calibration process of a polarimetric 
20 system involving the measurement of at least a reference sample in 
which 

• one illuminates the sample with a polarized incident light beam 
emitted by a polarisation state generator (PSG) containing a 
polarizer, said PSG modulating the light beam polarization, 

25 • said sample transmits or reflects a measurement beam, 

• one detects the measurement beam with an analysis section 
comprising a polarization state detector (PSD) containing an 
analyzer and detection means, and 

• one processes the electrical signals produced by the detection 

30 means with a processing unit. 

According to the invention, 

• Said PSG contains a first and a second liquid crystal elements 
LCj (j=1, 2) positioned after the polarizer, said LCj elements having 
retardations 5j between their ordinary and extraordinary axes and 

35 said extraordinary axes making angles 0j with respect to the 
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polarization direction defined by the linear polarizer so that by 
varying the retardation 5j of each LQ element for a fixed value of 
the 0j angle, when the LCj elements are nematic LCs, or by 
switching the orientation angle 0j when the LCj elements are 
5 ferroelectric LCs, one modulates the incident light beam 

polarization, the PSG having a modulation matrix (W) that is non 
singular, 

• Said PSD contains a third and a fourth liquid crystal elements 
LC’j (j=1. 2) positioned before the analyser, said LC’j elements 
10 being the same as the LCj elements of the PSG but positioned in 

the reverse order, so that by varying the retardation 6’j of each 
element for fixed values of 0’j angles LC’j when the LC’i are 
nematic LCs, or by switching the values of angles 0’j at fixed 6’j 
when the LC’j are ferroelectric LCs, one generates a detection 
15 matrix (A) for the analysis section, said matrix being non singular 

and so that for a given set of retardations (5j, 6’j) (j=1 . 2), or for a 
given set of orientation angles (0j, 0’j), one produces a measured 
quantity (Dn) and so that the processing unit produces the raw 
data matrix B = AMW, where (M) is the Mueller matrix of the 
20 sample. 

The processing unit produces after n=16 of such measurements and 
suitable data treatment: 

in a simplified (eiiipsometric) operation mode : the classical 
ellipsometric angles (T, A) as well as the overall transmission (or 
25 reflection) coefficient x characterizing the samples optically 

equivalent to DR, such as Isotropic non depolarizing surfaces 
measured in reflection. The measurement procedure includes a 
check of the validity of the description of the sample as a DR 
without the need of moving any part of the system, while with 
30 usual ellipsometers this can be checked only by rotating the 

analysis arm by 90“, 

in a complete (Mueller polarimetric) operation mode : the 
complete Mueller matrix (M) of any sample, with its sixteen 
coefficients, 

35 Said calibration processes thus comprise: 
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• for the simplified (ellipsometric) operation mode : 

- for ellipsometric measurements in transmission of samples 
assumed to be dichroTc retarders (DR), taking a complete measurement 
of a reference sample consisting of a DR defined by a Mueller matrix 

5 (Mo) with known parameters tq, To and Ao, said reference sample being 
propagation in air and (Mo) then being the identity matrix (lo), said 
measurement providing a reference raw data matrix Bo=AMoW, 

- for ellipsometric measurements in reflection of samples assumed 
to be dichroTc retarders (DR), taking a complete measurement of a 

10 reference sample consisting of a DR defined by a Mueller matrix (Mo) 
with known parameters (to, vi/o, Ao). said sample being a metallic mirror or 
a known sample for a system working in reflection mode (such as a NIST 
sample made of silicon covered by a known thickness of oxide), said 
measurement providing a reference raw data matrix Bo= AMoW, 

15 a) For the complete (Mueller polarimetnc) operation mode of a system 
working in transmission: 

- choosing a set of reference samples elements (p) comprising 
dichroTc retarders with approximately known Mueller matrices (Mp), 
defined by the parameters (xp.Tp.Ap) one of these elements being 

20 the identity matrix (lo) describing propagation in air, 

- for each of the reference samples (p), taking a complete 
measurement of said sample, set at an orientation angle 9p, by 
modulating the incoming light polarization and analyzing the 
outcoming light polarization, constructing the matrix (AR(- 

25 0p)MpR(0p)W) using the processing unit, this matrix being a 

product of the detection matrix (A), the Mueller matrix (R(- 
0p)MpR(0p)) of said element p set at the angle 0p, with R(0) a matrix 
describing a rotation by an angle 0 about the z axis and the 
modulation matrix (W), 

30 - calculating the product (AloW)'\AR(-0p)MpR(0p)W) for each 

reference sample p in order to obtain an experimental matrix (Cp) 
and determining Mp, or, more precisely, the values of its 
parameters xp, Tp and Ap, independently of the angles 0p through 
the eigenvalues of Cp, which are identical to those of Mp. This 
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allows a very accurate characterization of each sample in situ, 
during the calibration itself, 

- constructing a matrix (Ktot(9p)) equal to ^ (Hp(0p/Hp(0p)) 

P 

where the matrix Hp(0p) is defined as Hp(0p)[X] = R(-0p)MpR(0p)X- 
5 XCp where (X) is any real 4x4 matrix, 

- determining the eigenvalues Xi (i = 1 to 16) of the (Ktot(0p)) 
matrix in order to extract the modulation matrix (W) that verifies 
Ktot(W)=0, the p reference samples being chosen so that one and 
only one eigenvalue A,j vanishes when the angles (0p) used in the 

10 calculation of Ktot(0p) are set equal to their actual values during the 

calibration measurements, while the other eigenvalues A,j, being 
sorted in decreasing order of value, verify 7='k^5lX^ < 1 and the 
ratio Z is maximised, 

- This is equivalent to determine the modulation matrix W, 

15 (together will all the angles 0p), as the unique solution of the set of 

matrix equations: 

Mp(0p)X - XCp (8) 

20 - determining the detection matrix (A) by constructing the 

product (AloW) (yr\ 

According to various embodiments, the present invention also concerns 
the characteristics below, considered individually or in all their technical 
possible combinations. 

25 • a set of reference samples comprises 

- a linear polarizer set at 0i=O° orientation, 

- a linear polarizer set at 02 = 90° +/- 5° orientation, 

-a retardation plate with a retardation 5=110°+/-30° set at 
03=3O°+/-5°, 

30 • for spectroscopic applications, the retardation plate is an 

achromatic quarterwave plate. 

b) For the complete (Mueller polarimetric) operation mode of a system 
working in reflection: 

- choosing a set of reference samples comprising a linear 

35 polarizer, defined by its Mueller matrix Mpoi, and a first DR1 and a 
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second DR2 dichroYc retarders, said DRj having Mueller matrices 
(Mi), with i= (1, 2) respectively, with approximately known values of 
the parameters t \, 'Fi, Ai, 

- with each of the following sequence of elements, taking a 
measurement by modulating the incoming light polarization and 
analyzing the outcoming light polarization, the origin of the 
azimuthal angles (0=0) being taken in the plane of incidence, 

♦ DRi alone, set at 0=0, yielding a measured matrix 

Bi=AMiW 

♦ DR 2 alone, set at 0=0, yielding a measured matrix 

B2=AM2W 

♦ DRi, set at 0=0, and preceded by the polarizer set at an 
orientation angle 0i, yielding a measured matrix Bpi=AMiR(- 

01) MpoiR(0i)W, where R(0) is a matrix describing a rotation by 
an angle 0 about the z axis 

♦ DRi, set at 0=0, and followed by the polarizer, set at an 
orientation angle 02, yielding the measured matrix Bp 2 =AR(- 

02) MpoiR(02)MiW, 

- Calculating the products Ci=B 2 '^Bi and C 2 =BiB 2 '^ and then the 
matrices Ni=M 2 ’^Mi and N 2 =MiM 2 '^ through their eigenvalues, 
which are the same as those of Ci and C 2 . 

Ni and N 2 have actually the form of the Mueller matrices of DR, which 
are oriented, by definition, at 0=0. 

Calculating the products Cpi=B2*^Bpi=W^NiR(-0i)MpoiR(0i)W and 

Cp2=Bp2B2'''=AR(-02)MpoiR(02)N2A-\ 

Constructing a 16x16 real matrix Ki(0i) as Ki(0i)=H/Hi+ H pi(0i)^ 
Hpi(0i), where, for any real 4x4 real matrix X, Hi[X] and Hpi(0i)[X] 
are defined as Hi[X]=NiX-XCi and Hpi(0i)[X]=NiMpoi(0i)X-XCpi. 
Determining the modulation matrix W and the orientation angle 0i 
by requiring that Ki(0i) has one vanishing eigenvalue, and W is 
the vector associated with this vanishing eigenvalue. 

Constructing a 16x16 real matrix K2(02) as K2(02)=H2^H2+ Hp2(02)^ 
Hp2(02), where, for any real 4x4 real matrix X, H 2 [X] and Hp2(02)[X] 
are defined as H 2 [X]=C 2 X-XN 2 and Hp2(02)[X]=Cp2 X - X Mpoi(02)N2 




Determining the analysis matrix A and the orientation angle 02 by 
requiring that K2(02) has one vanishing eigenvalue, and A is the 
vector associated with this vanishing eigenvalue, 

reference samples are then chosen according to the following 
criteria: 

o the 16x16 real symmetrical matrices Ki(0i) and K2(02) will 
only have one vanishing eigenvalue, if and only if the angles 0i 
and 02 used for their evaluation are equal to the azimuthal 
angles of the polarizers during the calibration measurements, 
o The next eigenvalues are as large as possible, or, more 
precisely, the ratios Z=A,15/>.1 of the smallest nonvanishing 
eigenvalues (^15) over the largest (Xi) eigenvalues of Ki and K2 
are as large as possible. 

According to various embodiments, the present invention also concerns 
the characteristics below, considered individually or in all their technical 
possible combinations. 

• A set of reference samples is 

- a linear polarizer set at 0i=45®+/-5°, 

- a linear polarizer set at 02=-45°+/-5°, and 

- a couple of samples equivalent to a first DRi and a second 
DR2 dichroTc retarders, both oriented at 0=0 with respect to the 
incidence plane, with Mueiler matrices Mi and M2 such that the 
products M2 '^Mi and M2'^Mi are the Mueller matrices of a DR 
with 4^=45°+/-30° and A=90®+/-10, 

• For spectroscopic applications, said reference samples comprise 
a metallic mirror, 

• For spectroscopic applications said reference samples comprise 
an achromatic quarter-wave plate, oriented with one axis in the 
incidence plane placed before or after a metallic mirror. 

The invention also regards two measurement processes providing, after 
suitable instrument calibration 

• in the ellipsometric mode, the parameters (x,4^,A) of a 
sample assumed to be a dichroTc retarder (DR), 

• in the compiete polarimetric mode, the Mueller matrix (M) of 
any sample 




Both measurement processes involve, in all cases: 

• emitting an incident light beam linearly polarised along a 
direction of polarisation (i), 

• modulating the incident beam polarization, 

• sending the modulated incident beam to the sample, and 
returning a measurement beam, 

• collecting the measurement beam through a polarisation 
analysis section, 

• detecting the measurement beam after the polarization analysis 
section and producing electrical signals forming the raw data 
matrix B = AMW, 

• transmitting the electrical signals to a processing unit. 

According to the invention, 

• modulating the incident beam polarization by means of two 
liquid crystal elements LCj 0=1, 2) by varying either the angular 
orientations 0j, of the extraordinary axes with respect to the 
polarization direction (i) of the linear polarizer when the liquid 
crystals (LCs) comprise ferroelectric LCs (FLCs), or the 
retardations 6j at fixed orientations when the LCs comprise 
nematic LCs (NLCs), 

• producing measured quantities (□„) by means of an analysis 
section comprising two liquid crystal elements LC’j 0=1, 2) by 
varying the retardation 5’j of each element for fixed values of 0’j 
angles when the LCs are NLCs, or the values of the orientation 

angles 0’j for fixed values of retardation 5’j 0=1. 2) when the LCs 
are FLCs, 

The raw data B are then processed as follows: 

• in the eltipsometric mode : 

- calculating the matrix C = Bo’^ B = Mo'^ M W, where 
Bo=AMqW is the raw data matrix taken with the calibration sample. 
The eigenvalues of C are the same as those of Mq'^ M, which has 
the same form as the Mueller matrix of a DR. As a result, two of 
these eigenvalues (>,ri and A,r 2 ) are positive real, while the other 
two (^ci and Xc 2 ) are complex conjugates. 




deducing the ellipsometric parameters (x, T, A) of the studied 
sample from these eigenvalues and the known parameters xo, 'i'o 
and Ao of the reference sample according to ; 
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checking the validity of the description of the sample as a DR 
from the following relationship 



I'^cil (7) 

which must be obeyed by the eigenvalues of any Mueller matrix 
describing a DR, as it can be directly seen from eq. (3). 

• in the complete polarimetric mode : 

- the Mueller matrix M of any sample is calculated from the raw 
data matrix B as M=A'^BW\ 

To facilitate further the description of the invention, the following 
drawings are provided in which: 

Figure 1 is a schematic view of a polarimetric system operated in 
transmission according to the invention. 

Figure 2 is a schematic view of a spectroscopic ellipsometric 
system operated in reflection according to the invention. 

Figure 3 shows the experimental values of the diagonal block 
elements of the Mueller matrix of a high quality Babinet-Soleil 
compensator as a function of the setting x (mm) of the micrometric screw 
controiling the compensator retardation A (Fig. 3a), together with a plot of 
the values of A (deduced from the measured matrix elements according 
to eqs.(2)) versus the setting x (mm) of the micrometric screw (squares) 
and the corresponding linear regression (solid line) (Fig. 3b). 
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Figure 4 shows the off diagonal block elements of the same 
Mueller matrix versus the setting x (mm) of the micrometric screw 
controlling the compensator retardation A. 

These drawings are provided for illustrative purposes oniy and 
5 should not be used to unduly limit the scope of the invention. 

Figure 1 shows a polarimetric system according to an embodiment 
of the invention. It contains an excitation section 1 emitting a light beam 2 
and an analysis section 3. 

1 0 The excitation section 1 comprises a polarization state generator 4 

(PSG) through which passes the light beam 2. The polarization state 
generator 4 comprises a polarizer 5 that linearly polarizes the light beam 
2 along a polarization direction (i). 

First optical means 6 defines the geometry of the beam 2 at the 
15 sampie 7. 

The analysis section 3 comprises a polarization state detector 8 
(PSD) containing an analyzer 9 and detection means 10 for detecting the 
light beam 2. 

In a particular embodiment, the detection means 10 comprises a 
20 multipoint photosensitive detector that produces electrical signals sent to 
a processing unit 11. Said detection means 10 are adapted to 
polarimetric imaging and the multipoint photosensitive detector is 
advantageously a charge coupled detector (CCD). 

The polarimetric system may comprise as well a monochromator 
25 which is located in a first embodiment within the light source 12 that 
emits the light beam 2, before said beam enters the polarization state 
generator 4. In a second embodiment, the monochromator is located 
within the detection means 10, after the light beam exits the polarization 
state detector 8. 

30 According to the invention, the polarization state generator 4 of the 

polarimetric system comprises a first and a second liquid crystal 
elements 13 having birefringent axes, said liquid crystal elements 13 
being positioned after the polarizer 5. The polarization state detector 8 
also comprises a first and a second liquid crystal elements 14 having 
35 birefringent axes and positioned before the analyzer 9. The polarimetric 
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system comprises also control means for controlling said liquid crystal 
elements 13, 14. 

The present invention concerns Mueller polarimetric systems for 
analyzing samples represented by the sixteen coefficients of a Mueller 
5 matrix. The polarization state generator (PSG) 4 and the polarisation 
state detector (PSD) 8 comprise each a first and a second liquid crystal 
elements 13, 14 LCj (j= 1, 2) which may either be nematic liquid crystals 
(NLC) or ferroelectric liquid crystals (FLC). 

When NLCs are used, each NLCj element 13 of the PSG 4 

10 (respectively, for each NLCj element 14 of the PSD 8), has an 

extraordinary axis making a fixed angle 0j (resp. G’j) with respect to the 
direction of polarisation (i) and a variable retardation 5j (resp. 5’j) between 
its ordinary and extraordinary axes, which can be controlled 
electronically, said liquid crystal (NLCj) elements 14 being positioned in 
15 reverse order in the PSD 8 with respect to the NLCj elements 13 of the 
PSG 4. 

When FLCs are used, each FLCj element 13 of the PSG 4 

(respectively, for each FLCj 14 element of the PSD 8), has a constant 

retardation 6j (resp. 6’j) between its ordinary and extraordinary axes, and 
20 the angle 0j (resp. 0’j) between the extraordinary axis of the FLC and the 
direction of polarization (i) can be switched between two values 
separated by 45° by means of an electronic control device. Said liquid 
crystal (FLCj) elements 14 are positioned in reverse order in the PSD 8 
with respect to the FLCj elements 13 of the PSG 4. 

25 It is known that the application by control means of an appropriate 

voltage signal on each liquid crystal element (FLC or NLC) 13, 14 allows 
modulating the polarization of a light beam passing through said liquid 
crystal elements 13, 14. 

In mathematical terms, the liquid crystal elements LCj 13 of the 
30 polarization state generator 4 applies a polarization modulation such that 
the Stokes vector (S) of the light beam 2 at the exit of the polarization 
state generator 4 is given by: 
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where is the Mueller matrix of the LCj element (j= 1 , 2 ). When a set 
of four couples of retardations (61.62) or angles (01,02) are defined 
sequentially by the control means of the LCs, four linearly independent 
Stokes vectors are hence generated from an unpolarized light source. 

In a preferred embodiment, the liquid crystal (LC) elements 13 
according to the invention are nematic liquid crystal cells (NLC). Said 
liquid cells are particularly suitable to polarimetric imaging since their 
typical transmission range is currently between 400 nm and 1500 nm, 
and could be extended with new liquid crystal materials. 

With NLCs and in a preferred embodiment, the orientation angles 
0’j 0=1.2) are equal to 0j 0=1.2) and the retardations 6) 0=1.2) are equal 
to -6j 0=1.2) (modulo 27 t). Advantageously then, the couple of 
retardations (61, 62) takes sequentially the following values : (Ai,Ai), 
(Ai,A 2), (A2 ,Ai), {A 2 A 2 ), where Ai and A2 verify the formulae (315° + p 90°) 
and (135° + p 90°) respectively, where p is the same integer in both 
formulae and the angles 0i and 02 verify the formulae (e 27° + q 90°) and 
(e 72° + r 90°) respectively where e = + 1 has the same value in both 
equations while q and r are any integer, with tolerances on the angles 0i 
and on the retardations A equal to +-10° and +-20° respectively. This 
embodiment allows the simultaneous optimization of the condition 
numbers s(W) and s(A), both A and W matrices being non singular. 

The LC elements 13, 14 may also be ferroelectric liquid crystal 
modulators or any other suitable liquid crystal light modulating device. 

With FLCs and in a preferred embodiment, the extraordinary axes 
of the FLC make a coupie of azimuthal angles (0i,02) (resp. (0’i,02) ) with 
respect to the direction of the input polarizer 5 of the PSG 4 (resp. the 
output analyzer 9 of the PSD. and by means of a suitable electronic 
control device, these angles are set sequentially to (0i,02), (0i+45°,02), 
(01.02+45°), (0i+45°,02+45°). The azimuthal angles 0) are equal to 0j 
0=1,2) and the retardations 6) are equal to -6j 0=1.2) (modulo 2ti) for 
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simultaneous optimization of the condition numbers s(W) and s(A) of the 
modulation and analysis matrices. 

The retardations 5j are given by 6i=80°+-15° and 52=160°+-15°, while the 
orientation angles 0j are given by 0i = 67°+-10° and 02 = 160°+-40°. 

5 With FLCs and for spectroscopic applications (i.e. operation at 

variable wavelengths), as the values of the retardations (6i, 62 ) are not 
electrically controllable as for nematic liquid crystals, a birefringent plate 
is preferentially inserted between the two FLCs elements, leading to an 
overall optimization of A and W matrices in the whole transparency range 
10 of the FLCs, which is currently from 420 nm to 800 nm, and might be 
extended in the future with new materials. With typical values of FLC 
birefringence dispersion, the retardations (5i, 82) can be advantageously 
chosen equal to (90°, 180°) in the green part of the spectrum, while the 
birefringent plate can be chosen as a zero order quarter wave in the red 
15 part of the spectrum (633 nm), said plate being made of quartz. Hence, 
in a particular embodiment and for spectroscopic applications, the 
birefringent plate is a quartz plate and the PSG 4 comprises: 

♦ a linear polarizer, set by definition at an orientation angle 0=0, 

♦ a first ferroelectric liquid crystal, with a retardation 5i= 90°+-5° at 510 
20 nm, set at an orientation angle 01= -10°+-5°, 

♦ a quartz plate, providing a (true zero order) retardation 6q=90°+-5° at 
633 nm, set at an orientation angle 0q= 5°+-5°, 

♦ a second ferroelectric liquid crystal, with a retardation 62=180°+-15° at 
510 nm, set at an orientation angle 62= 71°+-10°. 

25 The invention regards as well conventional ellipsometry that is a 

special case of polarimetry for isotropic layers with smooth interfaces. 

Figure 2 shows a particular embodiment in which a spectroscopic 
ellipsometric system based on ferroelectric liquid crystals (FLCs) contains 
an excitation section 1 emitting a light beam 2, a sample holder 15 and 
30 an analysis section 3. 

The excitation section 1 comprises a polarization state generator 4 
(PSG) through which passes the light beam 2. The polarization state 
generator 4 comprises a polarizer 5 that linearly polarizes the light beam 
2 along a polarization direction (i). First optical means 6 focuses the 
35 beam 2 at the sample 7. 
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The incidence angle of the light beam 2 on the sample surface is 
defined as the angie at which the focused beam strikes the sample 
surface with respect to the normal to the surface 7. For example, a beam 
2 with normal incidence at the sample surface has an incidence angie of 
5 zero degree. The angle of incidence of the beam can be advantageousiy 
varied. The purpose of the focusing beam is to obtain a small spot on the 
sample 7, i.e. a compact spot with preferably dimensions inferior to a few 
tenths of mm^. This spot should provide a lateral resolution sufficient to 
map the sample surface. The light beam 2 emitted by the excitation 
10 section 1 is in the transparency spectral region of the FLC, which is 
currently from 420 to 800 nm, and might be extended further with new 
FLC materials. 

The beam 2 reflects off the sample surface and passes through 
the analysis section 3. In a more general case, the beam is scattered by 
15 the sample surface and passes through the analysis section 3. The 
analysis section 3 comprises an input optical (collimating) system 16, a 
polarisation state detector 8 (PSD) containing an analyzer 9 and 
detection means 10 for detecting the light beam 2. The detection means 
10 typicaily comprises a spectrometer coupled to several photodetectors, 
20 typically an array of CCD (charge coupled devices) that produces 
electricai signais. A processing unit 1 1 receives said eiectrical signals. 

According to the invention, the polarization state generator 4 of the 
polarimetric system comprises a first and a second ferroelectric iiquid 
crystai eiements 13 having birefringent axes, said iiquid crystai eiements 
25 13 being positioned after the polarizer 5, and a fixed retardation piate 17, 

between the two liquid crystal elements 13. The polarization state 
detector 8 aiso comprises a first and a second ferroeiectric liquid crystal 
elements 14 having birefringent axes and positioned before the analyzer 
9, and a fixed retardation piate 18 set between the FLCs. The 
30 polarimetric system comprises also control means for controlling said 
iiquid crystal elements 13, 14. 

The present invention can aiso be advantageously implemented 
for polarimetric imaging by a CCD camera. 

The polarimetric system and the polarimetric measurement 
35 process according to the invention have been the object of various 
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implementations whose following example demonstrates the quality of 
the results obtained. 

Figures 3 and 4 show the results obtained with a polarimetric 
system based on nematic liquid crystals, calibrated and operated in 
5 transmission, at 633 nm. The test sample was a high quality Babinet- 
Soleil compensator, which can be considered as dichroTc retarder (DR) 
with ^ « 45° and a retardation A which is a linear function of the setting 
of the compensator micrometric screw. The 16 elements of the 
compensator Mueller matrix (M) were measured for different settings of 
1 0 this screw. 

The experimental values for the diagonal block elements are 
shown on Fig. 3a, versus the screw settings (in mm) of the micrometric 
screw controlling the compensator retardation A; these values follow 
quite closely (to within 0.01) the behavior expected from eqs.(2) for 
15 4^=45°. In Fig. 3b is shown the variation of the dephasing A, as deduced 

from the values of the lower diagonal block elements, as a function of the 
screw setting: the standard deviation from a perfect linear fit is 0.13®, 
equivalent to A./2700, which is even better than the accuracy specified by 
the manufacturer of the compensator, thus confirming the performance 
20 of this polarimetric technique. 

Figure 4 shows the results obtained for the off-diagonal block elements 
of the Mueller matrix, again as a function of the micrometric screw setting 
(in mm). We recall that these elements are expected to vanish. Some of 
these elements, shown on Fig 4a, are independent of the compensator 
25 orientation, and they are always found to be smaller than 5.10'^ in 
absolute value. For the other ones (Fig 4b) this absolute value can reach 
1.5 10'^. This latter value might be due in part to an imperfect alignment 
of the components within the compensator itself. 




